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We study the mixing of the scalar glueball into the isosinglet mesons /o(1370), /o(1500), and 
/o(1710) to describe the two-body decays to pseudoscalars. We use an effective Hamiltonian and 
employ the two-angle mixing scheme for ri and rf. In this framework, we analyse existing data and 
look forward to new data into 77 and ?/ channels. For now, the /o(1710) has the largest glueball 
component and a sizable branching ratio into 7777 ', testable at BESIII. 
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I. INTRODUCTION 


Glueballs are arguably the most important uncon¬ 
firmed prediction of quantum chromodynamics (QCD) 
(see Refs. HH3] for reviews). Lattice QCD calculations 
predict the lowest-lying glueball around 1.6-1.7GeV, 
with quantum numbers J PC = 0 ++ . It is expected to 
mix with qq states in the same mass region, resulting 
in more 0 ++ states than naively expected from qq spec¬ 
troscopy. There seems to be some consensus that the 
states /o(1370), /o(1500), and /o(1710) are the relevant 
mass eigenstates made up from nn = (uu ± dd)/\/2 , 
ss, and the 0 ++ glueball G (see Tab. [I] for masses and 
widths). We note however that the lightest state /o(1370) 
is very poorly known, due to overlap with nearby states. 
To find the correct mixing structure one uses data on f 0 
both from production (seeking “glue-rich” channels) and 
decay; for this it is important to identify processes that 
are particularly sensitive to the glueball couplings. 

Gluons have a strong coupling to r] and rj 1 through the 
axial anomaly, which neatly describes the ratio T( J /ip -A 
?/7)/r(J/'0 -A 777) ~ 5 gH7]. For glueballs G - pro¬ 
duced for example in J/if -a G 7 - we expect T(G -A 
7777) > T(G -A ttk,KK) (expect the same from lattice 
QCD [HI S])- See Tab. [IT] for experimental decay rates 
of the potential glueball candidates /o(1370,1500,1710) 
into 7777 and 7777' from WA102 (soon to be updated with 
BESIII data ). 1 It is hard to see a coherent picture emerg¬ 
ing in the data. In particular, all /o(1370) data are incon¬ 
clusive/questionable to say the least. /o(1500) is proba¬ 
bly the best-studied state; only the 7777' channel is murky 
due to it being at the kinematic threshold. / o (1710) data 
come mainly from WA102 and desperately require confir¬ 
mation by BESIII. Our goal in this paper is to provide a 
framework to interpret upcoming data, which will hope¬ 
fully clarify the picture. 

The assignments of the three fo states to quark and 
glueball states are for the moment unclear, with different 
proposals in the literature [13]. Namely, Refs. na-EE] 
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suggested / o (1500) to be dominantly glueball in nature, 
which does not fare that well anymore [5]. Refs. [5] ITT] 
propose the mixing matrix 


/ /o(1370)\ 
/o(1500) 
\/o(1710)y 



0.78 

0.51 

—0.36X 

/ nn 

-0.54 

0.84 

-0.03 

I ss 

0.32 

0.18 

0.93 / 

V G 


(1) 


so / o (1710) is mostly a glueball (same qualitative picture 
found e.g. in Ref. m)- The decay fo -A 777 f is how¬ 
ever not discussed, because it is close to threshold for 
the /o(1500) - and hence subject to systematic errors - 
and inconsistently measured for the / o (1700). Since the 
glueball is expected to have a strong connection to the 77 - 
rf system, we will pay particular attention to these final 
states in this paper (see also Refs. mm)- 

A strong argument for the glueball nature of /o(1710) 
comes from the J/V’ decay rates [9j [ 21 ] | 22 ] 


T( J/y, -a 7/o(1710)) 
rm -4 7/0(1500)) 


10.5 ±6.5 


( 2 ) 


(where we used the tttt channel JS]), which is expected 
to be enhanced for glueballs due to J/ijj -A 'ygg -A 7 G. 
One also finds [23] 


T(J/V> -A7/o(1370)) 
T(J/V> -A 7/o(1710)) 


0.51 ±0.41, 


(3) 


which fits well into that picture. 

One more guideline for the identification of (predom¬ 
inantly) glueball states may be the total width, with a 
direct decay to quarks expected to be somewhat sup¬ 
pressed (by the OZI rule or large -N calculations). This 
argument is however difficult to implement at the level 
of two-meson decays, as the total cross section may very 
well be dominated by more complicated final states (in 
particular f 0 -A 47 t). Furthermore, some studies suggest 
rather broad glueballs [M], so we will not impose theo¬ 
retical conditions on the glueball width in this study. 
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State Mass [MeV] 

Width [MeV] 

a o (1450) 1474 ±19 

265 ± 13 

/o(1370) 1350 ±150 

350 ± 150 

/o(1500) 1505 ±6 

109 ±7 

/o(1710) 1722 ±6 

135 ±7 


g 547.86 ± 0.02 (1.31 ± 0.05) x 10“ 3 

g' 957.78 ±0.06 0.23 ±0.02 


TABLE I: Some relevant masses and widths [25]. 


Decay ratio 

Data 

Fit 1 

Fit 2 

r(a 0 (1450)^iCK) 

T(ao (1450)— tnrj) 

0.88 ±0.23 [25]* 

0 7c;+ 0 - 08 

u. /O-0.06 

o.8i!S;g? 

r(ao(1450)— ±Trr)') 
r(ao( 1450)— >7rrj) 

0.35 ±0.16 [13 

0.44 

0.44 

r(/ 0 (1370)^-7r7r) 
r(/o(1370)->XJC) 

12.5 ±12.5 [26]* 

11 + " 
1J --8 


r(/o(1500)—f7T7r) 
r(/o(1500)^iC7f) 

4.07 ±0.43 [25]* 

Q 00+0-94 

o.yy_ 0 65 

4-04±S;|| 

r(/o (1710)—>.7T7r) 
r(/ 0 (i7io)->xjc) 

0.4lt°;“ gg* 

0 46+ 0 04 
U ’^°-0.05 

0 20 +0 03 
U ’ ZU -0.03 

r(/o(1370)->r,Tj) 
r(/o(1370)—>-7T-71-) 

f) no+0.04 

U.UO_ 0 .03 

00 CO 

q q 
d 6 
+ 1 
CO 

0 

0 

0 04 +0,08 

U.Uf4_ 0#04 

r(/ n ( 1500 )— ►Tflrp 

r(/o (1500)—17T7T) 

0.145 ±0.027 [23* 

n 14+0.02 
u,1 ^-0.02 

0-14t°;° 3 

r(/ 0 (i7io)^ij7j) 
r(/o(1710)—17T7T) 

1 1 7+0-61 
± - ± * -0.48 

0 44 +0 10 
u -^-0.08 

0.96±°f 2 

r(/o(1370)-> W ) 

r(/ 0 (1370)^7fjf) 

0.35 ±0.30 [23 [28]* 

0.36l®j 6 

0-39 ±ttr 

r(/ 0 (i500)-.^?j) 

r(/o(1500)^iC7f) 

0.59 ±0.12 [25] 

0-56t“;« 

0.58l°;?° 

r(/o(i7io)-»7j»j) 

r(/ 0 (i7io)^iCif) 

0.48 ±0.15 [2K2Z]* 

0.20l“;°3 

f) 1 Q+0.03 

u.iy_ 0 .03 

r(/ n (i500)-W) 

r(/o (1500) yr/v) 

0.38 ±0.16 [23* 

0.53±g4? 

0 47+0- 14 
-0.12 

r(/ 0 (i7io)^^') 

r(/ 0 (i7io)-> w ) 

< 0.08 |27] 

0 64+ 0,22 

U.D+-0.22 

r» 09+0.20 

U.OZ-o is 

Br(/o(1500) — > 7771 ) 

0.349 ±0.023 [23* 

0 36+ 0 08 
VJ.OO_0.07 

n qc+0.07 

U.OO-0.06 

Br(/o(1710) -4 KK) 0.36 ±0.12 [22* 

rj q'7+0.06 
u -°‘ -0.06 

0 44 +0 07 
u *^-0.07 

Xmin/d.O.f. 


4 . 9/3 

6.4/3 


TABLE II: Data without a reference have been obtained by 
combining referenced data. Data with an asterisk are fitted, 
everything else are predictions. 


II. EFFECTIVE HAMILTONIAN 


To calculate the f 0 decays, we employ the effective 
Hamiltonian [50H52] 


H = h ± tr [X F PP] +h 2 G tr [PP] + h 3 G tr [P] tr [P], 

(4) 


with the pseudoscalar SU(3)f nonet 


P = 




7T 


K~ 


Vo_ 

v'S 


±2 ' V6 
K° 


Vo 

V3 


and the scalar qq states 


X F = 




K+ 

K° 

V6 + 



(5) 


( 6 ) 


See Fig.[I]for a diagrammatic representation of these op¬ 
erators (see also Ref. [13]). (Note that Ref. [33] proposes 
an entirely different effective Hamiltonian of chiral per¬ 
turbation theory; for the time being, we prefer not to rely 
on chiral perturbation in the energy range considered - 
more on this later about chiral suppression.) 

Assuming isospin invariance, we will adopt the stan¬ 
dard mixing scheme 


a o (1450)\ 

/o(1370) 

/o(1500) 

,/o(1710)/ 




V 



\ 

1 / 



(7) 


where U is a real orthogonal 3x3 matrix describing the 
mixing of nh, ss, and G into the /o states, further abbre¬ 
viated by (Fi,F 2 ,F 3 ) = (/o(1370),/o(1500),/ o (1710)). 

With the Hamiltonian from Eq. ([4]) one can easily cal¬ 
culate all the /o —» 7 T 7 T, KK, 7777 decays El. as well as the 
/o ~± w' decays of interest to us (which has been done 
already in Ref. [31], but only using hi^). The question 
is now how to include SU (3) /-breaking effects in Eq. Q 
- which should exist due to m s m u ( j. Some remarks 
are in order: 


• The coupling hi (Fig. [l] a) describes the OZI- 
allowed decay that should dominate the qq decays. 
The creation of an ss from the vacuum is modified 
relative to uu by a factor r a ~ 1 HZ], which we 
include explicitly in the amplitudes of hi. 

• h 2 corresponds to the dominant glueball couplings 
to the nonet (Fig. 0b). A flavour-democratic cou¬ 
pling of glueballs to quarks would lead to the ratios 

M(G -> 7T7T ) : M(G KK) : M(G -4 7777 ) 

:M(G -4 7 / 77 ') :M(G -a rpj') = 1 : 1 : 1 : 1 : 0 , 


and hence decay rates 3 : 4 : 1 : 1 : 0 modulo 
phase space 3J- Some of the relevant diagrams, 
however, contain the subprocess gg —» qq (dif¬ 
ferently attached gluons in Fig. [l] b) and it has 
been argued that this implies a chiral suppression 
cr(gg -4 qq) oc m q [32], 5B]- How this suppres¬ 
sion propagates to the hadronized G —► PP rate 
is unfortunately unknown. A proportionality of 
M{G —> PP) to current-quark masses is highly 
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FIG. 1: Diagrammatic representation of the effective couplings hi (a), h 2 (b), and h 3 (c). 


improbable for several reasons: (i) even accepting 
a Feynman-diagranr picture of these strong interac¬ 
tions, the gluons may land on different quark lines , 2 
(ii) even accepting the chiral suppression, we know 
that chiral symmetry is broken by confinement it¬ 
self, leading to constituent quark masses, (iii) al¬ 
ternate calculations based on perturbative QCD 
(equally out of their domain of application as chiral 
perturbation) suggest instead the decay constants 
fir, Ik , etc., as suppression Em mu, and (iv) com¬ 
pletely different arguments would use the overlap of 
wave functions at the origin as a factor governing 
the flavour nonuniversality . 3 

(Note also that the chiral suppression argument of 
the coupling of quarks to two gluons does not apply 
to the case of pseudoscalars through the quantum 
anomaly.) We thus leave the extent of the suppres¬ 
sion (if any) free for now, and introduce different 
couplings for the PP states to take the chiral de¬ 
pendence into account, following again Ref. HZ!: 

h 2 Gtr [PP] -¥ G [h 2 , n ((tt °) 2 + 27T+7T-) 

±2 h 2 , K (K°K° + K+K~) ( 8 ) 

+ ^2,o(?7o) 2 + ^2,8(?78) 2 ] ■ 

• h 3 (Fig. 0 c) describes the Gijorj 0 coupling that 
can be associated with the anomaly [301 [32 and 
is of particular interest to our study of 7777' final 
states. (In Ref. |39] this is taken one step fur¬ 
ther and promoted to a glueball-glueball-glueball 
coupling followed by a mixing of 77 and 77 ' with 
the 0 h glueball 03] ■) The full coupling is % D 
(/i2,o + 3/i3)G(77o) 2 , but we will fix h 2 fi instead of 
fitting it, so h 3 is indeed a free parameter for us. 

• Couplings tr [Xp] tr [PP], tr [X F P\ tr [P], and 
tr [X F ] tr [P] tr [P] would correspond to OZI- 
suppressed diagrams (see Ref. 32) and are ne¬ 
glected here as usual. 

• The free mixing matrix U of Eq. 0 describes in 
itself SU (3) / violation coming from the nondiago- 


2 Meaning cr(gg —» qq) oc m q is in any case only one possible 
subprocess. 

3 In Ref. 1381 . a coupling of the form MP 2 G is motivated within 

the Witten-Sakai-Sugimoto model for the singlet pseudoscalar 
go and subsequently extended to all P. 


nal mass matrix of the qq and G states. We will 
not assume a structure in this mass matrix (see 
e.g. Refs. 02j 03]) but rather determine U directly 
from the fo decays, keeping the mass eigenvalues 
and widths fixed to their PDG values (Tab. [i]) . 


A. 77—77' mixing 

Most studies of f Q decays assume the 77-77' system 
(Tab. [I| to be an orthogonal transformation of the 
SU(3)f states, 

fv\ _ (cos Op -sin 9 P \ (i l8 \ . . 

[rj'J (sin dp cos dp ) ( 770 ) ’ ^ ' 

with flavour eigenstates 770 = (uu + dd + ss)/a/ 3, r] S = 
(uu + dd — 2ss)/v / 6, and mixing angle dp = —11.4° [25] 
(see also Ref. 03]). 

It has been shown that this approach is a little too 
naive, because the large hierarchy <C together 
with renormalization effects, complicate matters pToTFfT] . 
In effect, one should rather use a general (invertible) 
transformation matrix, which has four real parameters, 

(v\ _ 1 (F s cos9 8 —Po sin0 O \ / r] 8 \ , 

\ pi'J F \Fs sin d 8 F 0 cos 9 0 ) yio) ' 

and is called the two-mixing-angle scheme . 4 (It reduces 
to the one-angle scheme in the limit F 0 = F 8 = P, 
9q = 9 8 = 9p.) Here, F = 92.2 MeV is a normalization 
constant, and a fit gives 08 ] 

P 8 /P = 1.26 ± 0.04 , 9 8 = (-21.2 ±1.6)°, (11) 
F 0 /F= 1.17 ±0.03, 60 = (—9.2 ± 1.7)°. (12) 

(Other methods of determining these parameters give 
similar values; see e.g. Refs 03HSZI-) 

Since this two-angle scheme is rather successful and 
stable in describing 77-77' data, we will keep the above 
parameters fixed to the best-fit values. For do (1450), we 


4 Using the quark basis g q ^ s instead of the flavor basis gs,0 in the 
two-angle scheme leads to two mixing angles 6 q , s that are acci¬ 
dentally close to each other, allowing (at current precision) for 
a one-angle description 6 q = Q s . Since this is obviously a basis- 
dependent effect, we will keep the two-angle formalism here. 
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then immediately arrive at the ratios 


r(ao(1450) —> 7 n/) 
r(ao(1450) —> nrf) 


( 2 F 8 cos 6* 8 + y^o sin 0 O \ p v > 
y 2F S sin 0 8 — V2F 0 COS 00 ) Pv 
0.44, 

(13) 


and 


r(a o (1450) -7 KK) 
r(ao(1450) —> nif) 


= / V^r a F o F 8 cos(0 o -0 8 )/.F 
V V 2 F 0 cos 0o — 2F$ sin 0 8 
— 0.84 , 



perfectly compatible with experimental results for r a ~ 1 
(Tab. |I5 • 

The /o decays are more involved due to the mixing 
matrix U, which we will determine by a numerical fit in 
the next section. The formulae for the decay rates are 
collected in the Appendix. 


B. Fit 


Having defined our Hamiltonian, we can perform our 
own fit to the decay data in Tab. [h[ using a simple y 2 
function. Note that we will keep the masses fixed (Tab. |Tj) 
and do not impose a structure on the /o mass matrix 
in flavour basis (given by t/ T diag (Mp 1 , Mp 2 , Mp 3 )U). 
Since there is unfortunately not enough reliable data to 
obtain a statistically significant result, our best-fit values 
have to be taken with a grain of salt. 

In order to calculate the rate /o(1500) —>• 77 ?/, which 
is at threshold because M v ± M v > = 1505.6 MeV > Mp 2 
(see Tab. [I]), we integrate the partial widths with a Breit- 
Wigner distribution: 

r {F 2 -> 7777 ') = \M(F 2 -> W)I 2 
T(F 2 -> 7777 ) \M(F 2 -± 7777 ) 1 2 

^ ImXmJ 2 d Ep(E,M v ,M v ,)/[(E 2 - M 2 ,) 2 + Mj, I*. ] 

X C;^ d E P (E, M v , M V )/[{E 2 - M 2 ,) 2 ± M^I%] ' 

(15) 

Here, 

, , J\E 2 - (Mi + M 2 ) 2 ]\E 2 - (Mi - M 2 ) 2 1 

p(E,M u M 2 ) = ^-— 

zh/ 

(16) 


is the momentum of a daughter particle in a two-body 
decay, where the initial particle has “mass” E. Since 
the F 2 — >• 7777' decay relies on a slightly off-shell F 2 , we 
integrated the 7777' and 7777 decay rates with the highly 
peaked Breit-Wigner distribution over [Mp 2 —Tp 2 ,Mp 2 + 
Tp 2 ], which reduces to [M v + M^,Mp 2 +Tp 2 ] for 7777' due 
to kinematics. Numerically, one finds 


T(F 2 -)■ 7777') I M ( F 2 ->■ 7777OI 2 

T(F 2 —>■ 7777) ■ \ A 4 ( F 2 —>• 7777 )| 2 ’ 


fairly stable against small changes in the domain of inte¬ 
gration. 

For the decay ratios of J/rj) we make the simple as¬ 
sumption that the decay is dominated by J/ip —» 77777 —> 
G 7 , leading to the simple formula 

T(J/ip ->■ 7 Fj) = \U i3 \ 2 p Fz . 

r(J /^^ 7 Fj) \U j 3 \ 2 p F / [ > 

but we will not include these ratios in our fits. See 
e.g. Refs. JS31E3] for dedicated studies. 


1. Fit 1: Flavour blind 


Taking for simplicity a flavour-blind G coupling h 2 = 
02 ,it = 02 ,k = 02 ,o = 02,8 an d the two-angle 77-77' values 
from above we can fit 0 i, 2,3 and the mixing matrix U, 
parametrized via three angles 9 t j (see Eq. (A2)). We 
find 5 


r a = 0.94^Q;g 5 , 0i = (698 ± 30) MeV, (19) 
02/0! = l-12±g;gt, 0 3 /0! = 1.5±8;1, (20) 

for the coupling constants and 

0 23 = 0.67 ±0.03, 0i 2 = 3.56 ± 0.04, 0 13 = 0.38^ ^ . 

( 21 ) 


for the mixing angles, resulting in the mixing structure 

(F 1 \ /-0.85 -0.38 0.37 \ (nn\ 

\f 2 = 0.53 -0.63 0.57 ss (22) 

\f 3 ) \ 0.02 0.68 0.73 / \GJ 

/ -0.91 -0.18 0.37 \ /singlet\ 

= 0.07 0.82 0.57 octet , (23) 

\ 0.41 -0.55 0.73 J \ G ) 


meaning that the largest glueball component is in 
/o(1710). The best-fit decay widths and ratios are given 
in Tab. HI 

We also gave the mixing matrix for the singlet-octet- 
G basis in Eq. (23), which seems to be a better starting 
point than nfi-ss-G, since the mixing angles in this ba¬ 
sis are smaller. / o (1500) is dominantly a flavour octet, 
similar to Refs. mm- While we have made no consider¬ 
ations on the mass mixing matrix, and have determined 
the mixing matrix U from decay properties, it is an inter¬ 
esting consistency check to verify to which original mass 
matrix M = U T M^gU in the flavour basis nh, ss. G 
our fit corresponds: 


/ 1393.09 -46.93 51.41 \ 

M = -46.93 1584.15 129.72 MeV, (24) 

y 51.41 129.72 1599.76 / 


5 The lcr range for parameter x has been obtained by fixing all 
other parameters to their best-fit values and then solving \ 2 (x) = 
Xmin + 1 for x - 


(17) 
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or 


/ 1412.53 -74.43 116.87 \ 

-74.43 1564.71 -76.24 MeV, (25) 
\ 116.87 -76.24 1599.76 ) 

in the singlet-octet-G basis. It is reassuring to see 
that the S'f7(3)-breaking elements are indeed smaller 
than the conserving (diagonal) ones and that no fine- 
tuned cancellations occur. However, the “octet mass” 
1565 MeV comes out surprisingly large compared to the 
states /\*(1430) and ao(1450). 

The ratios for J/t/> decays into 7 Fj show the qualitative 
behavior we expected from our initial considerations, 

T(J/^ 7/o (1370)) 

r m 7/0(1710)) ■ ’ 

r(j/y> -a 7/0(1710)) _ 
r(J/v> ->• 7/0(1500)) ■ ■ 

We find, however, a very small /q( 1370) width into PP, 


T(/ o (1370) -> PP) ~ 2 MeV, 


(26) 


dominated by the 7 T 7 t final state. From Eq. (A3) we see 
that this small rate is due to an interference of the nn 
and G amplitudes, reducing the rate by an order of mag¬ 
nitude. While we expect this rate to be subdominant to 
/o(1370) — > 4-7T [25], and the full width is known only very 
imprecisely, this might still be the most worrying result 
of our fit. (In Ref. [55. it is argued that the 47 t channel 
is subdominant to the 27 t channel.) 


2. Fit 2: Chiral suppression 

Since we expect a chiral suppression in the G -A PP 
decays, we perform a second fit with the ansatz 

*2,tt = *2 , *2,if = (fn/fn) 2 h2 , (27) 

*2.0 = ( Fo/F) 2 h 2 , h 2 ,8 = ( F s /F) 2 h 2 , (28) 

(motivated partly by Refs. |T71 lo7] l. which gives 

r a = 0.98/q q 4 , h x = (910 ± 40) MeV, 

(29) 

h 2 /h\ = —0.59 ± 0.02 , h 3 /h x = - 1.0 ±0.2, (30) 


for the coupling constants and 

023 = 0.60 ±0.02, 0i2 = 0.57 ±0.04, 0 13 = 0.51^2 > 

(31) 


for the mixing angles, resulting in the mixing structure 



0.74 

0.47 

0.49 

- 0.68 

0.55 

0.49 

-0.04 

-0.69 

0.72 

0.87 

0.04 

0.49 

-0.24 

-0.84 

0.49 

-0.43 

0.54 

0.72 



( singlet \ 
octet , 
G ) 


(32) 


(33) 


with only minor quantitative differences to the flavour- 
blind fit. The mass matrix in the nn-ss-G basis is given 

by 


and 


1421.53 

-48.29 

-61.19 

-48.29 

1574.06 

-143.85 

-61.19 \ 
-143.85 
1581.41 ) 

| MeV, 

(34) 

1426.85 

-55.81 

-133.02 

-55.81 

1568.74 

82.13 

-133.02 \ 
82.13 
1581.41 ) 

| MeV, 

(35) 


in the singlet-octet-G basis. In addition, we obtain the 
partial width 


T(/o(1370) -A PP) = 4 MeV, (36) 

and the ratios 

r(j/^^ 7 /o(l370)) 

T(J/V> -4 7 /o(1710)) ■ ’ 

7/o(1710)) _ 

T(J/V> -A 7/o(1500)) 

The same comments as above apply regarding the width 
of /o(1370). 


III. DISCUSSION AND CONCLUSION 

With the present state of data, the above fits can obvi¬ 
ously not convincingly decide on the main issues of glue- 
ball spectroscopy, namely the composition of the /o states 
and their decay properties. We consider the present 
attempt mainly as a preparation for the interpretation 
of future data. Nevertheless, in both fits we find that 
/o(1710) has the largest admixture of the 0 ++ glueball, 
whereas /o(1500) comes out close to the octet flavour 
structure. The very small rate of /o(1370) -A PP is 
particularly concerning, even though / o (1370) -A 47 t is 
expected to make up most of the width. The coupling of 
the glueball to r/o plays an important role in both fits and 
is numerically large, making rj and r{ final states crucial 
testing grounds. We certainly expect / 0 (1710) -A 7777' to 
be visible in BESIII, which will hopefully clarify some of 
the issues. It may also be surprising to some that we 
obtain a better fit with flavour-blind glueball couplings. 

We stress that the use of the two-angle scheme for 77- 
rf mixing is important when discussing the /o -A 7777 (') 
decay modes, as the usually employed one-angle scheme 
is inconsistent. 
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Appendix A: Decay rates and ratios 

The 3x3 orthogonal mixing matrix U that describes 
the mixing of the /o states and the qq and glueball states 
is defined by 


with Cij — cos Qij and Sij = sin 6 ij. For the decay rates, 
we define p 2 ,x = h, 2 ,x/hi and p 3 = h 3 /h\, leading to 


, 1 3 /i? 

r (FW — 


(V%Uii + 2p2,nUi3^ ^2 • (A3) 


/ /o(1370)\ [FA (nn 

/o(1500) = F 2 = U ss 

\/o(i7io)y \f 3 \g 


Here and below, p\ — \Px\i where p x denotes the three- 

(Al) 

momentum of X in the two-body decay F t —t XY. In¬ 
troducing 


and for the numerical fit parametrized by three mixing 
angles 0 ± 2 , 623 , and 9 13 using the common PDG notation 


Ai = V2Uu + Fr a Ui2 + 6 / 32 , 8^3 1 


Bi = Un - \/2r a Ui2 , 


/ C12C13 

S12C13 

Sl 3 \ 



I —C23S12 — S23S13C12 

C23C12 — S23S13S12 

S23C13 I 

Ci = 2 

VZUn + r a Ui2 + 3 (p 2 ,o + 3p 3 )Ui 3 

V S23S12 — C23S13C12 

—S23C12 — C23S13S12 

C23C13 / 




(A2) we find the ratios 


(A4) 

(A5) 

(A6) 


T(Fj -» 7T7T) 

r(Fi -A KK) 

r(F ->■ 7777 ) 
T(Fi 7T7r) 

r (Fj -a irn') 
T(Fi -4 ryq ) 


3 / V%Un + 2p2 }7r Ui3 \ p„ 

2 yr a Un + \/2Ui2 + 2\/2p2,K'Ui3 ) Pk 

1 ( F 1 2 3 4 5 6 7 8 9 [F$Ai cos 2 9 0 - AF 0 F 8 Bi cos 9 0 sin 6 8 + F 8 2 C,; sin 2 0 8 ] \ Pr} 
27 \ F 2 F 2 cos 2 (d 0 - 0 8 ) [V2U xl + 2 p 2 ,nU. l3 \ ) V*' 

1 / F^Aj sin 2d 0 + 4 F 0 F 8 B X cos(g 0 + g 8 ) - F 2 Q sin 2g 8 \ 2 

2 V Fq Ai cos 2 9q — 4FoF 8 -Bj cos 9 0 sin 0 8 + F^Ci sin 2 0 8 ) p v 


(A7) 

(AS) 
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